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Abstract 

We assume that the Stefan problem with undercoohng has a classical so- 
lution until the moment of contact of free boundaries and the free boundaries 
have continuous velocities until the moment of contact. Under these assump- 
tions, we construct a smooth approximation of the global solution of the Stefan 
problem with undercooling, which, until the contact, gives the classical solu- 
tion mentioned above and, after the contact, gives a solution which is the 
solution of the heat equation. 

1 Introduction 

In this paper, we study the confluence of free boundaries in the Stefan problem 
with kinetic undercooling in the spherical symmetric case. 

In the preset paper, we study the problem in the domain Q = Q x [0,ti], 
where J7 = [i?i,i?2] is an spherical layer in the spherical coordinates. We 
assume that the domain Q is divided into three layers ^f2{'t) and ^~{t) as 
follows: 

nt{t) = [Ri,n{t)], n-{t) = [n{t),r2{t)], n+{t) = [r2{t),R2], 

where i = 1,2, are the free boundaries of phases "+" and "— ". We 

assume that the phase "+" occupies the layers (7^2 (*) and the phase "— " 
occupies the layer Q~{t). 
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We shall construct a smooth approximation of solutions of the Stefan 
problem with undercooling under the assumption that the motion of the free 
boundary is the motion of the front of a nonlinear wave and the confluence 
of free boundaries is interaction of solitary nonlinear waves. The possibility 
of this interpretation is given by the models of phase field [2] proposed by 
G. Caginalp. In fact, the choice of the method for approximating the limit 
Stefan problems with undercooling is unessential for us, because we do not 
prove that the approximations thus constructed are close to the corresponding 
solutions of the phase field system. 

For example, we could use the definition of the generalized solution of the 
limit problems including the order function (the nonlinear wave) [lOj . 

Here our considerations are based on the following simple fact. Suppose 
that there are two families of solutions (exact and approximate solutions) of 
some problem depending on a small parameter e. Suppose that both these 
families have the properties that permit passing to the limit as e — > in the 
weak sense. Suppose also that the family of approximate solutions satisfies 
a problem with a right-hand side small as e ^ in the weak sense. Then 
the weak limits of both families are solutions of the same limit problem and 
if the latter has a unique solution, then the difference between the exact and 
approximate solutions tends to as e — >■ at least in the weak sense. 

We recall that the phase field system has the form 

L9 = -^, eLu - ^^—^ - 6 = 0, (1) 



where 



d I d f 



The function = 9{r,t,e) has the meaning of the temperature, and the 
function u = u{r,t,£), which is called the order function, determined the 
phase state of the medium: u ~ — 1 corresponds to the phase "— " in the layer 
Q~{t), and u ~ 1 corresponds to the phase "+" in the layers r2^2(^)- 

Passing to the limit as e ^ in ([1]), we obtain the Stefan problem with 
undercooling. 

This passage to the limit is possible, for example, in the case where the 
corresponding limit problems have classical solutions. In this case, the weak 
limits as e — > of solutions ([T]) give these solutions [21 (Sj [H] . 

We denote 

Oir,t) = ^. (2) 
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Then the phase field system ([T|) we can write in the form 



da d'^a _ du . 
u-u^ a 

eLu = h-. (4) 

e r 

The smooth approximations of solutions of the Stefan problem with ki- 
netic undercooling constructed in this paper are approximate (in the above 
sense) solutions of the phase field system and they admit a weak passage to 
the limit as e ^ 0. In this case, we obtain the limit problems (and their so- 
lutions) describing the process of confluence of free boundaries. To construct 
these approximations, we use the assumption that the classical sharp fronted 
solution of the Stefan problem with kinetic undercooling exists until the con- 
fluence of free boundaries begins. This is a natural assumption in our paper 
(otherwise, it is not clear the confluence of what is considered) and can be 
provec0. 

We repeat that, under this assumption, we construct an approximation 
of solution of the limit problem, which is smooth for e > and uniformly 
bounded in e ^ 0. As is well known, in a similar problem about the propa- 
gation of shock waves, the existence of such an approximation distinguishes a 
unique solution. 

By t* £ [0,ti) we denote the instant of confluence of free boundaries. 
Then, for any t ^ t* — 6 for all 6 > 0, we see that the limit problems have 
solutions. The asymptotic solution of system (l3]),(|4]) has the form 



= a 



(r.t) + (S+(r,t) - S-(r,t)) u,, (^^^L^"j |'!jt|(*)\ _ 



V — f(t) 
+ LOO I ^ I + e 



£ 



r + f{t) r — r{t) 



Here uji{z) — > 0, 1 as z — > =Fcxd, uj^^\z) € S(Kz) for A; > 0, Lp{t) is a smooth 
function, ujq{z) = tanh(z), and uj{t,zi,Z2) G C°°([0,r]; ^(M^)). By S(M") 
we denote the Schwartz space of smooth rapidly decreasing functions. If the 
initial data for ([T]) has the form ([SD, (0) at f = 0, then, for t ^ t* — 5, we have 
the estimate 

\\u - ur ; C(0, T; l2(m1)|| + ||a - aT]C'{Q)\\ ^ ce^, /x ^ 3/2, 



"'^A. M. Meirmanov. A private communication. 
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where {a, u) is a solution of system ([3]), ([3]) (see [HIS]). Here Q = x [0, t* — 5), 
and the constant c is independent of e. 

The main obstacle to the construction of solutions of the form ([5]), ([6]), 
which could be used to describe the confluence of free boundaries, is the fact 
that, instead of an ordinary differential equation whose solution is the function 
^0(2) [21 E], in the case of confluence of free boundaries, we must deal with a 
partial differential equation for which the explicit form of the exact solution 
is unknown. 

In the present paper, we use the technique of the weak asymptotics method 
[H E], which allows us to avoid this problem. Let us explain several basic 
points. 

Definition 1. A family of functions /(r, t, e) G L^{Q) integrable with respect 
to r for all t G [0,ti] and for e > admits the estimate OT>'{e'^) if, for any 
test function C(r) G C^, we have the estimate 



Generalizing ([7|), we shall say that the family of distributions f{r,t,e) 
depending on t and e as on parameters admits the estimate Ox)'{£'^) if, for 
any test function C{f)i we have the estimate 



Example 1. Let u>{z) £ S(M^), Jj^i = 1, and let xq £ [^1,^2]- Then we 

have 



Example 2. Suppose that tJi(z) G C°°, (wj)' G S(M^), lim^^_oo ^i(^) = 0, 
and lim2_>+oo "-^il-z) = 1, i = 1,2. Then we have 

(a) 





where the constant C^^^^ depends on ti and the test function C{f)- 



(/(r,t,e),C(x)) =0(e^), O^t^h. 





(b) 




Bi( ] H{x-xi) + B2{ H{x - X2) + Ov'ie) 
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where Ax = xi - X2, Bi{p) £ C7°°, B'-ip) e 8{R^), Bi{p) + Bi{p) = 1, 
Bi{+oo) = 1, and Bi{-oo) = 0. 

Example 3. The preceding relations easily imply the formula 

H{x - xi)H{x -X2) = B [^f^ H{x - xi) + (1 - B)H{x - xa) + Ov'{e), 

where B[p) £ C°° is the function Bi from the preceding example. 
Example 4. The following relation is a corollary of Definition 1: 

for all m £ Z4. and q > 0. Generally speaking, this is not true for the 
derivatives w.r.t. t. 

The relations in Examples 2(a), 2(b) and 4 are obvious, but, for the 
reader's convenience, we prove some of these formulas in Section 5. 

Here we only note that, in view of items (a) and (b), one can say that 
the product of approximations of the Heaviside functions or of the Heaviside 
functions themselves is a linear combination with accuracy up to small terms 
(~ C'x)'(e)). The general nonlinear functions of linear combinations of approx- 
imations of the Heaviside functions can be linearized similarly. Precisely this 
property underlies the constructive study of the interaction between nonlinear 
waves with localized fast variations. 

Definition 2. A pair of smooth functions {u, a) is a weak asymptotic solution 
of the phase field system ([3]),(j4]) if, for any test functions C{^)i(,{^) ^ C'^($7), 
the following relations hold: 

/ {rut + at)r\dr+ [ ^.(r^O.dr = 0(e), (8) 
Jn Jn 

e / UrUtidr + 2e / riil^dr - - / ul{r'^^)rdr (9) 
Jn Jn 2 

" 7 I (t " T + i) + X u{rai\dr = O(e^), p e (0, 1/2). 

The left-hand side of Eq. ([9]) is obtained by multiplying equation ^ by 
and integrating by parts. The reminders 0(e"), a = l,p, in the right-hand 
sides of dS]) and Q must be bounded locally in t, i.e., for t G [0,ii], we have 

max |C'(e)| ^ Ctj^e, Ct^ = const. 
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This construction was introduced and analyzed in [6]. 

The fact that a noninteger exponent appears in the right-hand side of ^ 
is not directly related to the technique of the weak asymptotics method, see 
Examples 1-4. The source of this noninteger exponent is the nonsmoothness 
of the function 9 (the temperature), which appears at the instant of confluence 
of the free boundaries. 

The paper is organized as follows. 

In Section 2, we explain the structure of ansatzes of approximations of the 
temperature and the order function. These ansatzes are constructed under the 
assumption the classical sharp fronted solution of the Stefan problem exists. 
At the beginning of Section 2, we formulate what we exactly need. 

Next, in Section 3, we substitute the constructed ansatzes in system ([3]), ([4]) 
and derive equations for the unknown functions contained in the ansatzes. 

The results of these calculations are summarized in Theorems 1-3 in Sec- 
tion 3. The final result of this paper is formulated in Theorem 4 in Section 3. 
It states that the assumption on the existence of the classical solution is 
sufficient for constructing formulas for the weak asymptotic solution of sys- 
tem (ISl),® in the sense of Definition [TJ Next, in Section 4, we analyze the 
constructed formulas and derive the following effects: 

(a) the weak asymptotic solution is smooth for t > t*, the absolute values 
of the free boundaries velocities are equal to each other at the contact 
moment; 

(b) the temperature has a negative jump at the instant and at the point of 
confluence of the free boundaries, and this jump is equal to 

[^11.=.*, = -^^ lim (r?oi+rio,). 

where r = rjo(t) other positions of the free boundaries till the interaction. 

In particular, it follows from (a) and (b) that the velocities of the free 
boundaries have jumps at the point of contact. 

These effects can also be discovered in the numerical analysis of the process 
of confluence of free boundaries. Some helpful technical results are given in 
Section 5. We note that the only example known to the authors, where the 
confluence of the free boundaries is studied, is given in |llj . 

In general, this paper turned out to be technically rather complicated and 
long in spite of the fact that the details in several justifications in Section 5 
were omitted. 

Although, as was shown above, the results obtained here are a necessary 
step in the study of the multidimensional problem, this paper shows that the 
following classification can be introduced: 
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(1) problems in which the confluence of free boundaries leads to disap- 
pearance of one of the phases; 

(2) problems in which the domain occupied by one of the phases changes 
its connectivity, but the number of phases remains the same. 

In this paper, we consider an example precisely from the first class of 
problems. As was mentioned above, we do not justify the asymptotics of the 
constructed solution. But this can be done based on our constructions. It is 
easy to see that the main role here is played by the (not proved) existence of 
the classical solution up to the moment of confiuence of the free boundaries. 
Our construction of the existnece under this assumption reduced the justifi- 
cation to estimating the soluton of the heat equation with the right-hand side 
fe admitting the estimate 



and with zero initial and boundary conditions. An analysis of the structure 
of this right-hand side shows that /g • as e — > is a linear combination of 
functions 6'{r — rt) and S{r — ri), i = 1,2 with coefficients depending on t, r, 
and these coefficients converge fast to zero as r ^ itcxD. 

Hence we can conclude that for r ^ ri the solution of this heat equation 
belongs to C°° for e > and admits the estimate 0(e'^). In the whole domain 
Q X [0,T], a rough analysis based on general theorems [13] shows that the 
solution belongs to W2^, 6 > 0, and has the estimate 0{e'^) in the norm of 
this space. In particular, the weak limit of the constructed weak asymptotic 
solution is equal to the exact global solution of the heat equation in the phase 
field system. 



2 Ansatz of the approximation of the solu- 
tion of the Stefan problem with undercool- 
ing 

In our case, the Stefan problem with kinetic undercooling has the form 



/iG (0,1/2) 



da d'^a 



r G [Ri,R2] 



r^hit), i = l,2. 



(10) 



dt dr'^ ' 




(11) 



Ida! 




(12) 



dr 



J r=fi{t) 
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Relations (|10p - (jl2p are supplemented with (Dirichlet or Neumann) boundary 
conditions for r = Ri, i = 1,2, and with a consistent initial condition. 
Obviously, problem (jl0p - (jl2p can be written as 

f\f2t^{r -f-i) -r\fxt^{r -fx), (13) 



dt dr"^ 
a 



r 



r=fi{t) 



We shall assume that the initial conditions are chosen so that the problem 
in question has the classical solution, r2(0) — fi(0) > 0, and there exists a 
t = t* £ [0,ti] such that fi{t* — 0) = f2{t* — 0). More precisely, we assume 
that 

(i) the limits limt^t*~o'rit, i = 1,2 exist; 
(ii) 

n {c^ {n-it)) u i^tit)) u {n+{t))) n c{n) 

(here D denotes the interior of the domain D), cf. [3j. 

As will be shown below, these assumptions permit describing the interac- 
tion (the confluence of free boundaries) constructively. 

We shall construct this approximation (ansatz) as a weak asymptotic so- 
lution of system (l3l),([l]). 

First, we introduce the ansatz of the order function in (ll0p -( fT2|) . It has 
the form 



1 + a;o ( P^-^] +CO0 ( P^-^ ] (14) 



ri-r\ „ r-r2 
- wo I 13 — - — I ujQ ( (3 — - — 

Here the unknowns are the functions (3 and = rj(t, e), i = 1, 2. In what 
follows, we write them more precisely. Now we only note that if /5 > 0, then, 
for ri < r2, the function u coincides up to 0{e^) with the sum of the first 
three terms in the right-hand side of ()14p , for ri = r2 , we have u = 1 + O-jj' (e) , 
and for ri > r2, we have u = 1 + 0{e^). These relations can easily be verified 
directly. We note that, under the above assumptions, we can continue the 
functions rj(t) to the interval [0,ti] preserving the smoothness and the sign 
of the derivatives. We choose such continuations and denote them by r.to(t). 
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Now, we can write the functions ri{t,£) and (3{t,e) more precisely. Namely, 
we set 

ri{t,e) = rio{t) +i;Q{t)rii{T,t), (15) 

where 

tpo{t) = r2o{t) - rio(t), T = ipo/e. (16) 

We note that 

r ^ CO, t < t*, 
T — > — oo, t > t*, 

as e ^ 0. Thus, the value of the variable r characterizes the process of 
confluence of free boundaries. Also, we introduce the functions 

i>{t,e)=r2{t,e)-ri{t,e), p{t) = iP{t,e)/e (17) 

We set 

/?(r) = l + /3i(r). (18) 

In addition, we assume that 

/3j")(r), ^^)(r)^0(r-i), r ^ oo, a = 0, 1. (19) 

< 6i < /5i(r) < 62, 61, 62 = const, 

rn ^ r:r^, m - r:^^ = 0{\t\'^), r ^ -00, r^:^ = const. 

Assumptions (jl9|) for the functions rn are the assumption that as e — > the 
functions ri{t,£) approximate continuous functions with a possible disconti- 
nuity (jump) of the derivatives at t = t* . 

More precisely, the functions rj(t, e) determined in (jlSp are approximations 
of the functions 

ri(t,0) = rio + ipor'Hi-ipo), 

where r^^ = limr^-oo ?"ii(''", t), ri{t,e) — ri{t,0) = 0{e), and H{z) is the Heav- 
iside function. 

Now we describe the ansatz of the approximation of the temperature 
a{r,t). Here we also must use the continuation procedure for matching the 
temperature for t < t* , when there are three subdomains of the domain Q, 
and the temperature for t > t* , when there is only one phase. For this, we 
first introduce a "model" of the temperature whose continuation is reduced to 
the problem of continuation of functions depending only on the time t. This 
model must have the same structure as the temperature a. Obviously, the 
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simplest function of this form is a function linear in x in the layers 0^"^^, {t) and 
quadratic in the layer Q~(t). 
Namely, we set 



da 

dr 
da 
dr 



= ±7f , 



r=ri±0 



r=f2±0 



The functions 7,^ = "yfit), j = 1, 2, depend on t, are continuous for t < t*, 
and have limits as t ^ t* — 0. Therefore, they can be continued to the interval 
[0,ti], ti > t*, so that the properties of being continuous and the signs are 
preserved. We shall use the previous notation for the continued functions. 
Now we can introduce the temperature model 

f = Ji{t){ri - r)H{ri - r)+j+{t){r - r2)H{r - ra) 

+ ^-(r,t)^Jl^^^}^l:ilAH{r - n)H{T2 - r) (20) 

W 

+ ^^r, - T)H{r - rs) + /(r, t), 

where 7~(r, t), 7(r, t), and I{r,t) are linear functions of r, ^~\^^^. = 7j"(i), 
7\r=ri = 7i(*)) and I\r=ri = (" l)*"'"'^(»^jnt + >^2)- In more detail, the functions 
7(r, t) and 7~(r, t) can be written in the form 



7 

/ = 



2 

71 +72 
2 

ki - k2 



- (r - r*) 

— (r — r*) 



71 - 72 

ki + /C2 



2 ' ' 
We have the relations 

ri{t,e) =fiit) + 0{e 

This implies that, for t <t*, 



ri + r2 



ki = ririt + X2, k2 = r2r2t + X2- 



1,2. 



dT 
dr 

df' 
dr 



dl_ 
dr 



(21) 



(22) 



5{r - ri) 



r=ri 

da 
dr 



5{r-ri) + Ov'{s). 
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We rewrite the product 



as follows 



H{r — ri)H{r2 — r) 



H{r - n)H{r2 - r) = B{p)[H{r - n) - H{r - rs)] + Ov'{e). (23) 
H{ri - r)H{r - rs) = (1 - B{p))[H{r - rs) - H{r - n)] + Ov'{e). 

Here B{p) G C°°, B{p) ^ 1, p ^ oo, B{p) ^ 0, p ^ -cx), = 0(1^1"^), 
IpI — > oo for any > 0, a > 0. 

In what follows, we write B{t) instead of B{p). The difference is that B{p) 
is an unknonwn function [p is unknown), while B{t) is a known function. 

But we do not use this remark and, by definition, we set 



/ + 



-7^" (n - + — (''I - - r2)B{p)- 



-(ri - r)(r - (1 - B{p)) 



H[r — ri) + 



7^(ri - r) ^(''i ~ '^)('' ~ '''2)-B(p) + 



(ri - r)(r - r2) (1 - B{p)) 



H{r — r2). 



Now we can write the expression d'^T /dr"^ uniformly in -0 = r2 — ri. We 
have 

= 1 [-i^{t)B{p) + 71 (t) (1 - B{p))] H{r - n)+ 



[72-(t)i?(p) +72(t) (1 - B[p))\ H(r - r2) + 



(24) 



btii) + 7r(t)i?(p) - 71 (i) (1 - B{p))\ 5{r - n) + 
btit) + 72it)B{p) - Ut) (1 - Bip))] 6{r - r2), 

We set = —71 and 7^ = 72 and see that the coefficients for 6{x — ri) in 
(l24l) take the form 



[71" + 7i+] B{T)6{r - n) + [72 + 72+] B{T)6{r - n] 



(25) 



We note that the replacement of B(p) by B{t) implies a correction of order 
C'x)'(e). Now we note that the inequalities 



[H{r - ri) - H{r - r2)] /i> > 0, [H{r - r2) - H{r - n)] /tp < 0, 
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hold for ip ^ 0. Also, we note that dT/dr is a smooth function. So, by (j22p 
and the Stefan condition (|12p we have 

7r + 7i+ = r?o(t)r'io(t), 72" + 72"" = -rloitywit) 

for t < t* and assume that the continuations of the functions contained in 
these relations are chosen so that these relations hold for t £ [0,ti]. 

The assumption that a(r, t) is the classical solution of problem (jl0p - (jl2p 
implies that 

a-f = qeC{[0,t*);C\n)) (26) 

and the limit q exists as t — > t* — and for x £ fl. Therefore, we can continue 
the function q to the domain 17 x [0, ti], where ti > t* and the properties of 
smoothness are preserved, and moreover, 

Q\r=n=0, i = l,2, t<t*. 

We shall construct the global temperature 9 in the form 

a = e{r)f + q + q, (27) 

where q is the desired function, e(r) G Cq°{[Ri, R2]), and e = 1 for r G 

[n(0),f2(0)]. 



3 Construction of the weak asymptotic so- 
lution 

We first consider the heat equation ([3]). Using the formulas of the weak 
asymptotics method and taking (fT4|) into account, we obtain 



9u ^ 2 



+ 2TTri 



ru{2 - B^,) + 



P' 



00 



-r2,(2-i?oo) + M.5f 



^2 ^00 

where the estimate Ox>'{£) is uniform in r, ipQ = r2o — rio. 



5{r-ri) (28) 
S{r-r2) + Ov'{e), 



B, 



00 



ujo{z)ujo{-r] - z)dz, B 



00 



zcjo(^)t^o(-?/ - z) dz, 77 = p/3. 



We assume that function q and its continuation are chosen so that the 
function q satisfies the boundary conditions of the original problem. Then the 
boundary conditions for the function q are zero. 
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We substitute the function a determined by relation (p7j) and expression 
dMD for du/dt into Eq. i^. We denote 

1- — - — 
dt dr"^ 

According to the definition [2] with accuracy C'x)'(e), we obtain 
dii 



La + r— = LeT + Lq + Ai6ir - n) + A2S{r - rs), 
at 



where 

q = q + Q, 



i = l,2. 



P2 -00 

We let T tend to oo (i.e., for t < t*) and, in view of (10) and (18), obtain 

L{ef + q) = f|r2(5(r — f2) — r|r2(5(r — fi). 

Since eT + q ^ a as t ^ oo, we see that q ^ as r ^ oo. The total 
equation for the function q = q + q has the form 

L{q + e{x)I) = F{x,t,T) 

d'^ ( ^{ti — r){r — 



Qj.2 
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B{T)[H{r -n)- H{r -r2)] (29) 

- 1^ (^^^Il^lR^lA^j (1 _ i?(r))[i/(r - rs) - /7(r - n)], r G 

where -F(r, t, r) is a piecewise continuous function containing of terms that 
are uniformly bounded in e on O x [0, fi]. 

Equation was derived under the assumption that 



By qi, q* we denote the terms in q such that 
2(7i+ + 72" 



Lqi 



-B{T)[H{r-r2)-H{r-n)], 



(30) 



(31) 



Lq2 = ^^^^^^^(1 - B)[H{r - ra) - H{r - n)], 

Lgi* = 2(r - r*)(7+ - j^)i;~^B{T)[H{r - rs) - /7(r - n)], 
Lg* = -2(r - r*)(7- - 72+)V^-2^(T)[i7(r - rs) - /7(r - n)]. 
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With accuracy up to functions smooth in $7, we can calculate them us- 
ing the fundamental solution of the heat equation (in what follows, we shall 
consider only the first equation), 

2v27r Jo V* - a 
-exp{-(.-e)V(t-a)}^^^^_ 

It is clear that the functions qi are uniformly bounded. It is also clear that 
Qi i C^(17), but qi E Ci'2(J7 X [0,ti]) \ ({r = n} U {r - r^}). Therefore, to 
justify Eq. ([29j) . we must verify that no (5- functions arise in calculating the 
derivatives d'^qijdr^, d'^qljdr^ and dqi/dt, dq*/dt. 

For this, we note that, for any test function ^{r) up to functions smooth 
in 0,, we have the relation 

{QiXir)) = ll\i{a)B J%-'f{C,t-a)da]dC, (33) 



In 

\2 



dr is the solution 



where /(e,t- a) = /iRi C(r-)e-(''-«) / ^2T:{t - a) 

of the heat equation ft — f^^ = CiO the point t — a. It is clear that 
/(^,t) G C°°(0) for all t. Calculating the derivative dqi/dr, we obtain 

(| ^ , = , C) = 2 ^* Tf {a)B j'J ij-^h (e, t-a)dad^, 

where fi{(,,t) is a solution of the equation ft — f^^ = C"{0- It is clear that 
the relation 

'||,C(r))^EG^C(r.) + o(.). 

cannot hold for any coefficients Gi. Similarly to (j33|) . with accuracy up to 
smooth functions, we can write the functions q^, i = 1,2. They have the same 
properties as cji with the only additional condition 

qi\r=r* = 0. 

This easily follows from the explicit formula of the type of (31) and the fact 
that, for r = r* , the integral with respect to ^ is an integral of an odd function 
over a symmetric interval. Therefore, condition (j32p is necessary for deriving 
Eq. ()29p . This implies that the function ^ + eJ is uniformly bounded in 
X [0, ti] and belongs to C^'^ for r / r^, z = 1, 2. 
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So if the function q + el satisfies Eq. (j29p with zero initial conditions 
and the boundary conditions that follow from the fact that the function 6 = 
eT + q + q must satisfy the boundary conditions of the original problem, while 
the functions T and q are known (more precisely, they will be determined 
after the functions (pi are found), then we have the following assertion. 



Theorem 1. Suppose that the function a is determined by relation {27), and 
the function q is a solution of Eq. i29\) with zero initial condition and zero 
boundary conditions. Suppose that relation \30\) holds. 

Then the pair of functions a and ii is a weak asymptotic solution of the 
heat equation in the phase field system, i.e., relation ([^j holds. 



Remark 1. Expression (j30p may seem to be strange, because we did not 
take into account the well-known fact that the Dirac functions are linearly 
independent. But this fact is taken into account in the framework of the weak 
asymptotics method (see Lemma 2, Section 5). Here we do not want to mix 
the substitution of the ansatz into the equation and the analysis of the results 
of this substitution (see Section 4). 



Now we consider the Allen-Cahn equation (|4]) 
weak asymptotics of the expression 



We must calculate the 



dcf 



du 
dr 



eLu 



[ii — 



where the function a is determined by the relation (j27p . We obtain 

j:=Vl5{x-^i) + V^5{x-^2) (34) 
+ V^5'{x + Vi5'{x - ^2) + Ov\e), 

where V-, i,j = 1,2, are linear combinations of several convolutions. Their 
expressions will be given below. Here we note the following. We have 



Uj.{r'^Ur)rC,dr -\ — j u,. (-u — n^) Qdr 



(35) 



e 

2 ./rol 



where 



Fiu) 



n4 



n2 



1 

^ + 4- 



We denote 



1 



^(^' ^) = n {1 + ^0(2) + t^o(-z -if)- ujq{z)u:q{-z - r/)} . 



(36) 
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We have the following estimates: 

Q{z,r]) = uJo{z) + fi{z,r])e~'^'', 

where 



-oo, 

r/ ^ oo. 



(37) 
(38) 



\f{z, ri)\dz < const, 



1/1(2;, T])\dz < const. 



These relations readily follow from the explicit form of the function ujo{z), 



UJo{z) 



+ e 



In fact, relations (j37l) and (j38|l express the above-described properties of the 
ansatz u (jl4p in different terms. 

Using the technique of the weak asymptotics method (see Lemma [7] in 
Section 5), we can write 



/ r^rUtCdr = V rfnt [ ^'^{^', - n' ) dz 

J ^',{^',-Vl'^)dz + Ov'{e) 



(39) 



/32 



Similarly, we obtain (see Lemmas [5] and [71 Section 5): 



rUrCrQdr = —ri {el + q + < 



\r=r\ 



^r,{z,r])C{ri)dz 



(40) 



+ r2 (el + g + 



\r=r2 



^vi^, v)C{r2) dz + O(e^), fx G (0, 1/2). 



Thus, adding expressions (j35|) . (j39j) and (|40|) . we see that the coefficients 
of the (5-functions in formula (1341) have the form 



Vi 

where 



ruBn + 



TJZ 



+ i-iy^'u {el + q)l^^. Cn - nC, (41) 



i = l,2, Bn = j ^'M'z-%)dz, 

Bl, = j [z{^'z - -{z + - f^;) dz, (42) 

Cn = J in', - n'^) dz, q = q + q 
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and 



1 



c = - I (n'^Ydz (43) 



Similarly, we obtain (see Lemma [6l Section 5): 

V,^ = rl (^I3C + ^Dj , Vi = rl (^PC + ] , (44) 



where 



b = -^ I F{n)dz (45) 



Thus, to obtain P = (3{r]) we have the equation (obviously, this is a neces- 
sary condition for the relation T = O-jy^e) to hold, see (134]) ): 

/?2 = £ (46) 

According to (f43ll and (f45ll {D is positive), the right-hand side of the relation 
dM]) is positive. 

So we have proved the following assertion. 

Theorem 2. Suppose that the assumptions of Theorem 1 and are satis- 
fied and 

2 

J2 ViS{x - fi) + Vi6{x - ^2) = OD'{e). (47) 
1=1 

Then the pair of functions a, u is a weak asymptotic solution of the phase 
field system ^ in the sense of Definition\^ 

Thus, under the assumption that the classical solution of the phase field 
system exists (see (i) and (ii) above), relations (j29]) . (pOj) . (06]), and (JIT]) are 
sufficient conditions for constructing a weak asymptotic solution of system 
([3]), (Jl]). In what follows, we prove that the relations mentioned above are 
equations for determining the functions [3 = (3{t) and rji(r), i = 1,2. We 
present an algorithm for solving these equations. 



4 An analysis of the relations obtained 

We begin with relation (|30p and, for a while, forget everything said about the 
notion of linear independence. 
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Then, for this relation to hold, it suffices to have the two relations 

BV^t+li\=^^^ ^ = 1>2. (48) 

In view of (j22|) . we have +7j~ = [^] Taking this into account and 

adding relations ()48p and Lemma [21 we obtain 



1 

~ 2 
B{t) 



lo{z)i^q{-z -rj)[pr- 2^z 



de 

dr 



+ 



r=fi(t) 



89 
dr 



r=f2{t) 



dz 



(49) 



Here we used the relation 

ipt = r2t - nt = ip'oPr + ^"0(^^21 - rn)t, 

which follows from the definition of the function p, see (|17p . 

We set the last term to be zero, since we show below that r2i — rn = 
Odrl"^). In view of Lemma 3 in Section 5, in this case we have (r2i — 
i^ii)t4'o = C(£)- Moreover, in view of the Stefan conditions and the choice of 
the continuation of the functions and roj(t), we have 



+7i +72 



89 
8r 



+ 



89 
8r 



r=fi(t) 

We transform the left-hand side of relation (1491) : 



r=r2{t) 



T def , / 
I = Vo 



ioQ{z)u]Q{-z -v) { Pt- ) dz 



Pr + ^ 



uJo{z)uJo{-z - r])^zdz 



djo{z)uJo{-z -'n)\ pT + -02^-^ -V) + -^V ] dz 



Next, we change the variables in the last integral 



/ {-z-r])i;jo{z)uJo{-z-r])dz 
and note that 



zuiio{—z — 7])uio{z) dz 



1 8 
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Finally, we get 



1 

J_d_ 

2pd^ 



{zujoiz)uJo{-z -r]) + zujoi-z - r])ojo{z)} dz 
Coo{z)uJo{— z — T]) dz . 



Since the function ujq{z) is odd, we see that the expression in braces in the 
first integral in the right-hand side is 

d 

-Z— (1 - LJq{z)u}q{z + T])) . 

Hence, integrating by parts, we obtain 



2/32 
1 dr] 
2^8^ 



(1 - uJo{z)uJo{z + r/)) dz 
LOoiz + ri)uJo{z) dz . 



Or, finally. 



where 



Pt + 



1 d f\ 



2dT \p 



Biv) 



(1 — uJo{z + r])ujQ{z)) dz = 2r]tanhr]. 



So, in view of (j49p . we obtain the equation for p: 

1 



d_ 

dr 



P+^^B{r,) 



Or 



d / rj + r] tanh r] 



After integration, we obtain 

r/(l + tanh??) = 2/3 



2B{t). 



2B(t). 



B(a) da. 



(50) 



(51) 



(52) 



Since B{t) — > 1 as r ^ 00, we have r]T~^ — > 1 as r ^ 00. 
Since B{t) = 0(|t|~^) for any as r ^ —00, the integral in the right- 
hand of (j52p converges and hence the left-hand side tends to the limit 

lim r]{l + tanhr]) = lim 2(3 / B{a)da = 0. 
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Moreover, in view of the inequality B ^ 0, we have ^ for r G M^, which 
impUes that 

r/ — > as T ^ — oo. 



Here we took into account the inequahty (3 > 0, which follows from (j46p . 

Substituting /? = {CD~^{r]))^^'^ into (j52p . we obtain the following equation 
for the function r]: 

7?(1 +tanhr?) = ((7I)"^)^/2 /" (53) 

Jo 



Since the functions contained in (j53p are monotone and the limits exist as 
r — > ±00, this equation is obviously solvable. 
Next, we have 



/3(r) = 4 / £^ (3~ = const, r ^ -00. 

Moreover, as readily follows from the exponential rate of convergence to the 
limit of the functions C/D, 1 — B, and 1 + tanh?7 as r — > ±00, the derivatives 
satisfy the estimate 

^ = 0{\r\-^), \t\^oo. 

Q^a VI I II 

Now we calculate the limits of the expressions 

Ji = B{-if + -ir)-A„ i = l,2, 
as T — > ±00. As r ^ 00, we have 

p^oc (p^t), B^l, Ai^i-iy+'rlr',^. 
Therefore, we have 

lim Jj = 

in view of the Stefan conditions (jl2p . and Jj = 0{t~^), i = 1,2. 
As r — > — 00, we proved that r] ^ and hence 

ri-r2 = 0{e), r ^ -00. (54) 

Theorem 3. Relation ^5^^ is a sufficient condition for relation ^3^) to hold. 
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Proof. In view of the corollary of Lemma [2l Section 5, the estimate for the 
coefficients of the (5-functions in (j30p as r ^ oo and ()54p are sufficient condi- 
tions for relation (|30p to follow from (09]). In turn, relation (|53p follows from 

(BSD. □ 

Now we analyze the relations V-^ = and (I47p . We use the explicit form 
of the function -F(-z) and the formula 

+ e ^ 

As r/ ^ oo, we have 

n{z,r]) =coo{z) + Oie~^'^) (56) 

and hence /? ^ 1 as r/ — > oo. 

From relation ()36p we easily obtain 

0;(z,r/) =0(e-l^l), |z|^oo, 



l^;(z,r?) = 0(e-l^l), |z|^oo, 
Cn = \ i^o{z){l - uJo{-z - 7])) dz > 0. 

Similarly to (|56p . using (j48p . we can verify the relations 

/3, = 0(e-2H), r/^oo, 

= l + 0(e-'l''l), r/^oo, (57) 
0^ = 4 + 0(6-2"), ry^oo, 
5^ = 0, r/ ^ oo. 
We rewrite the expressions for the coefficients V^^ in more detail: 

V,' = -rfntBn - ^-^^Bl, + n (el + q)^,^ Cn - nC, (58) 

Vi = -rir^tBn - ^^^^^ - (el + q)^,^ Cn - r^C. (59) 

We take the limit in formulas (I58p . (I59p as r ^ oo i.e. prior to interaction. 
So, in view of conditions (03]) . ([57|) we obtain 



r=rio(t) no(i)' 



r=r20(t) r2o(i)' 



21 



where = J ujQ{z)dz. In view of (jlip we obtain xi = X2 = C~^. 

Thus, in view of the corollary of Lemma [2l Section 5, for the relation 

- + V^6ix -ip) = Ov'{e) 

to hold, it is sufficient to have /r\ + V2 /r"^ = or, in more detail, 



{rit + r2t){Bn + C^)-'^-^Bl, 



+ 



r=ri 



r=r2 



/32 



(60) 



1 1 

ri r2 



^2Cn — C 



0. 



Thus, we have proved the following assertion. 

Theorem 4. The conditions of Theorem 1 and relations i53\} . |57p , i60\) are 

sufficient for the functions a, u to be a weak asymptotic solution of the phase 
field system Q). 

We consider relation ()60p . It follows from the above that Cn = C^{ri) 
decreases sufficiently fast as — > 00, in any case ||fy|CQ(??)| ^ const. This 
fact, the estimate 

q\r=ri - q\r=r2 = 0(|n " ^'2^), G (0, 1/2). 

proved in Lemma[3]in Section 5, relation (j54p implies that the obvious estimate 
I In - ra^Col ^ e>'const\r]^'Cn\ = 0{e^') 

for t ^ t*. 

For t ^t*, by Lemma El we have 



'r=rio 



where -01 = '"21 — ''iij Hence, by (jl7p . and a statement similar to Lemma [3l 
in Section 5, we have 



2 (riot + r2ot) 













+ (7+ 


(i 




r=rio r 


r=T-2o 


\rio 


^20/ 



Oie^"), t ^ t* 



We introduce a function V{t) G C°°, such that £ S{R^), V{-oo) = 0, 
V{oo) = 1. Then, in view of considerations similar to those used in Lemma 3, 
we can show that the following estimate hold: 



{Bn + Cn) V{T){riot + raot) + 
+ (^2Cq - C 













r2j 



0{e' 



+ - 

r=ri r 



r=r2 
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This implies that the left-hand side of is estimated as Ox>'{e'^) if 

{Bn + Cn) [{riot + raot) + nt + rat] - 2^5^ = 0. (61) 
From this relation we obtain 

|M.„+r.)) = -(n„.+r.„,) + 2^M_S_^ (62) 

From this equation we determine the extensions of the functions rn and r2i. 
We note that our argument results in an equation that does not contain the 
temperature, namely, the assumption that the classical solution exists till 
t = t* is sufficient for constructing a global smooth approximation of the 
solution. 

Let us calculate the function Bq = BQ{rj) in more detail. We have 

Bn{v) = I n',{n', -n'^)dz. 

In this integral, we make the change of variable z ~z — t]. Then we obtain 

o;^-(o;-o;), (o;-o;)^-o;. 

Hence we have 

Bn{7]) = J Q'.n'^dz. 

Finally, we obtain 

Bn = ll {n'M - ^',) + ^'z^'r,) dz = \j {n'^f dz. 

Thus, for finite r], the denominator {Bq + Cq) in the last term in the right- 
hand side of ()62p does not vanish. Moreover, using the explicit form of the 
functions ioo(z), B^, Bq, and Cq, we can verify that 

niPrBn /n/i„i-iv 



where ^ 1 is an arbitrary number. 
Therefore, we have 
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We note that 

V'o (l + ^[^(^21 - rn)]) = VoPr, (64) 

where p = 7](3~^ and the functions r] and (3 are determined by Eqs. (|53|) 
and dM]). 

The system of Eqs. ([631) and ([Mil allows one to find the functions rji, 
i = 1,2, and thus completely determines the functions contained in the ansatz 
of the weak asymptotic solution of system ([3l) , dl]) . 

It is easy to see that the solutions thus constructed satisfy conditions ([19]) . 

Next, using the explicit formulas for B^, Bq, and Cq, we can easily verify 
that B^\jj=Q = (because uJo{z) is odd) and (-BQ + Cn)|^=o Oi because cjo(-z) 
is even. This means that 



(rit + rat) 1^=0 = 0. (65) 

We note that the right-hand side of ()29p is a piecewise smooth function for 
{ipl ^ const > and smooth for r 7^ r^, i = 1, 2. As ip ^ {t ^ t*), the right- 
hand side of (f29l) becomes proportional to (5(r— Tq), where Tq = r*\^=Q, and the 
proportionality coefficient is negative. Thus, if r and t vary in the respective 
neighborhoods of the points Tq and t*, then the (5-function with a negative 
coefficient appears and disappears in the right-hand side of (j29p . which results 
in a negative soliton-like jump of the temperature in a neighborhood of the 
point t = t*, r = Tq. 

Calculation of the temperature jump. To prove the statement in the Intro- 
duction concerning the temperature jump, we must calculate the quantity 

[el + q + q]\t=t*,r=r*- 

Let us verify that g + g is a continuous function. This function is the sum 
of solutions of the heat equation with singular right-hand sides, but these 
singularities arise after the substitution of the continuous function eT (which, 
of course, is not a solution) into the left-hand side. It remains to prove that 
the solution differs from the function eT by a continuous function. For this, it 
suffices to verify that the right-hand sides, which arise after the substitution 
of eT into the equation, do not generate any singularities in the solution of 
the heat equation in addition to those contained in the function eT. 

Here we, in contrast to the preceding statements, use the fact that the 
singularities in the right-hand sides arise as the result of the substitution. 

By Qi, i = 1,2, we denote the terms in q+q corresponding to the right-hand 
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sides: 



f - {'Pi-x){x-(p2) 



B{T)[H{r-n)-H{r-r2)], 



h = - (^T^^^^^l^) (1 - Bir))mr - r.) - - n)], 

The other terms in q arc solutions of the heat equation with piecewise con- 
tinuous right-hand side and hence are continuous. 

The functions qi have a similar property if they are calculated up to 0(e). 
Indeed, we denote 

- (n - r)(r - r2) 



n = 7" 



and represent, for example, the function qi in the form 
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1 I'BirV.e)) rB^-^^^ifg,^,,,^ 



Omitting the number factor and integrating by parts in the integral over ^, 
we obtain 



-f 

Jo 

-f 

Jo 



B 



+ 



Vt-t 

B 

d 



■ e 4(t-t') 



r2 



e 72 - e" 



dt' 



dt' y/t^ ' 



The last term is a solution of the heat equation with the piecewise contin- 
uous right-hand side 

B^[H{r-n)-H{r-r2)], 

and hence it is continuous uniformly w.r.t e > 0. The other terms are also 
continous functions uniformly w.r.t e > 0. 

It is easy to sec that in the formula for qi there is no term containing the 
derivative because this term would be of order 0{e). 

Indeed 

dB^ 



—Il[H{r-n)-H{r-r2)] 
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since |n| ^ Ctp for r G [ri,r2] and e ^ip = p. The derivative B'^ decreases 
faster than any power of |/o|~^, < const. This imphes that 

Hence we have 

[q + q]\r=r*,t=t' = 0. 

Let us calculate the function el. We have 



It follows from (ISD) that 



Pr ^l, r ^ oo, 
Pr ^ 0, r ^ — oo 



and hence 



1 



[el]\r=r*, = -T ^ lim (r^ot + raoJ 



So, the jump of the temperature has the form 

[^11.=.., = -"-^ ^ hm (rfoi + rL). 

5 Technique of the weak asymptotics 
method 

First, we recall the definition of regularization of the generalized function. 

Definition 3. A family of functions /(x,e) smooth for e > and satisfying 
the condition 

w — lim/(2;, e) = f{x) 

e^O 

is called the regularization of the generalized function/(x). 

We note that, by definition, the last relation can be rewritten as 

lim(/(x,e),C(x)) = (/,C) 

for any test function ^(x) (from now on, ( , ) denotes the action of a generalized 
function on a test function). 
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Lemma 1. Let Tt = {x — ip{t) = 0}, x G M, where ip{t) is a smooth function, 
let co{z) £ S (§ is the Schwartz space), and let (5 = f3{t) > 0. Then the 
following relation holds for any test function C{x)- 

^(^ (/?^^) , C{x)) = ^AUiv) + 0{e), 

where = j^^uj{z)dz. 

Proof. The expression in the right-hand side can be written as 

i / Jf-^^)ax)dx=^-^ I .iz)dz + 0{e). 

Here we perform the change of variables z = P{x — (p)/e and apply the Taylor 
formula to the integrand at the point x = ip. By definition, the last integral 
is the action of the generalized function P~^A^6{x — ip) on the test function 

C- □ 

Proof of the formula in Example 1 in the Introduction. Let uj{z) £ 
C°°, J £ §(IR^), lim^^+oow(z) = 0, and lim^^_oo ^(2:) = 1- We verify that 

CO ( ) — H{x — xq) = Ojy/{e), Xq = const, (66) 



e 

where H is the Heaviside function. By Definition[Tl we consider the expression 
uj ^ — — H{x — Xq) C{x) dx (67) 

= eaxo) I [u;{z)-H{z)] dz + 0{e^). 

Here we performed the change of variables z = {x — x^j/e and applied the 
Taylor formula to the functions C,{x) at the point x = xq. In view of our 
assumptions and the properties of the function the integral in the right- 
hand side of the last relation converges and hence the right-hand side of (j67p 
is of order 0(e). We thus obtain estimate ([66]) . 

Proof of the formula in Example 2 in the Introduction. Let uJi{z) G 
C°°, (uJiY e S(Mi), lim,^+ocUJi{z) = 0, and lim^^_oo a;i(^) = 1, i = 1,2. We 
consider the integral 

jde{ f fx-xi\ fx-X2\^, 
J = wi W2 C{x) dx 



1 \ e / ye 

1 [ ) ~~P^ ' ' ^ C(y)'i2/ ) dx. 



r C{y)dy) 

J -00 J X 
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Integrating by parts, we obtain 



dx. 



We perform the change of variables z = {x — Xi)/e with z = 1 in the first 
integral and with i = 2 in the second integral and apply the Taylor formula to 
the function F{x) = J^^Civ) dy at the points x = Xi, i = 1,2, respectively. 
Then we calculate the first and the second integral, we have 

/ Ax\ / Ax\ 
J = Bi{ H{x -xi) + B2{ H{x - X2) + Ov'{e), Ax = xi - X2, 



Bi{p)= LOi{z)uJi{-z - p)dz, B2{p) = uj2{z)ui{z - p)dz. 

To calculate the linear combination of generalized functions up to Cx>'(^"); 
we must improve the classical definition of linear independence. This improve- 
ment plays the key role in the study of problems with interaction of nonlinear 
waves. 

Indeed, let 0i / ^2 be independent of x. We consider the expression 

gi6{x - 4>i) + g25{x - ^2) = a > 0, (68) 

where the functions gi are independent of e. Clearly, for the last relation to 
be satisfied, it suffices to have 

5i = 0(£«), i = l,2, 

or, with the properties of the functions gi taken into account, 

5i = 0, i = l,2. 

But, if we assume that the coefficients gi depend on the parameter e, then 
the above estimates do not work. Namely, let us consider the following specific 
case of this dependence: 

gi = (t>'}Si{/\<l)/e), n>l, z = l,2, (69) 

where the functions Si{a) decrease sufficiently fast as \a\ — > 00, and = 

(j)2 - (j)!. 
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Lemma 2 (Linear independence of generalized functions). Suppose that the 
estimate 

\aSi{a)\ ^ const, i = 1,2, —oo < a < +00 
holds. Then, for a = 1, expression ()68p implies the relation 

Si = -S2. (70) 

Proof. Using the Taylor formula in ([68]) and taking ([69|) into account, we 
obtain 

<A?5iC(<Ai) + cP^S^Cih) =riSiQ{ct>i) + </.^52C(<Ai) 

+ <A2^2(<A2-</'l)C'(01+/^02), 

where < < 1 Since the function aS2{(j) is uniformly bounded in cj G R^, 
we obtain 

52(A<^/e)(</.2 - </>i) = {-a52(a)}U=A0/, • e = 0{e). 
If Si = —S2 then expression (f68]) can be rewritten as 

{cp^Si + 0^S2)C('/'i) = (</>? - </'^)5i = 0{e). 
Thus we obtain the assertion of Lemma [2j □ 
Corollary 1. Suppose that 

|cS'i((7)| < const, i = 1,2, a >0 

and the functions 4>i = (j)i{t,e) are continuous and uniformly continuous and 
satisfy the condition that Ac/) > for t > and e > and At/) = 0{£) for 
t < 0. Then, for a = 1, expression (66) implies (64). 

Proof. Following the above argument, we obtain 

SlCiM + ^2C(02) = (<Al - cl>2)Sl + <P^S2ih - <Pl)C'i<Pl + f^H). 

In view of our assumptions, 

</'^52(<A2 - <Al) = 0(e), (</)? - </)^)5i = 0(e) 

uniformly in which implies the desired assertion. □ 
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Lemma 3. Suppose that f{t) G C^, f{to) = 0, and /'(to) 7^ 0. Suppose also 
that g{t, t) locally uniformly in t satisfies the condition 

\Tg{t,T)\ ^ const, \Tg'i-{t,T)\ ^ const, — oo < r < oo, 

and g{tQ,T) = 0. Then the inequality 

\g{tJ{t)/e)\^eCi, 

where Cf = const, holds in any interval ^ t ^ t that does not contain zeros 
of the function f{t) except for t^. 

Proof. The fraction is locally bounded in t. The fraction '^f^l^^'' is also 
locally bounded. We have 

g{tj{t)/e) fit) t-to 



g{t,f{t)/e)=e 



(t - to) e fit) • 



According to the assumptions of the lemma, the last factor in the right-hand 
side is bounded on the interval under study. The product of the first and 
second factors (without e) is bounded in view of the properties of the function 
g{t,T). □ 

Corollary 2. Suppose that the assumptions of Lemma [3] are satisfied for 
^ T < oo (— oo < T ^ 0). Then the assertion of Lemma^ holds on any 
half-interval {to,t\, which does not contain zeros of the function f{t), and 
sign£= sign/(t), t e {to,t]. 

The proof of Corollary [2] is obvious. 

Lemma 4. The following inequality holds: 

q\r=ri - q\r=r2 = 0(^1 " r2\^)- 

Proof. We consider only one of the the functions q from (j32p . namely, the 
functions determined by relation (j33p . and apply the method developed in 
We consider the difference of expressions (|32p omitting the number 
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factors: 



\r=ri H\r=r2 

I ('•l-fl' , ('-2-0% , 

Jo Jrh - T 

i p exp( j[exp IJ 



7^:4^ i^Vt - 1' 

{t-t'Y {r2-n)^' 



(r2-ri)M {t-t'] 



^ exp(- ~ )fexp ~ ^^^^"^ + ~ - 11 

Jti IpVt - t' 



where we choose fj, G (0, 1/2). 
Next, using the inequahty 



< const 



for a > 0, r £ [0,oo), and taking into account the fact that the integral 



(t _ t')-(l/2+M) rfi'^ (0,1/2), 







converges, we obtain the statement of the lemma for r2 ^ ri. If r2 ^ ri (this 
is true for t > t*), then we must change the exponents outside the square 
brackets in the integrands. □ 

Lemma 5. Suppose that u>{z) G C°° decreases faster than any power of \z\~^ 
as \z\ 00. Then 

e~^Lj((r - ri)/e)qj = qj5{r - ri) / io{z) dz + ©©/(e^). 



where qj is one of the functions defined by the relations [31\) . i = 1,2, /x G 
(0,1/2). 
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Proof. We consider one of the functions q, namely, the function defined by 
the relations 



1 



dt' 



2\/vr7o Wt-t 
The desired relation can be written as 



exp( 



0' 



C(r) 



r - n 



qdr = C,{ri)q\ 



uj{z)dz + Oie^") 



or, omitting the number factors, in the form 

■(r - ri)' 



C(rV 
dt' 



qdr 



(71) 



i^Vt- t' Jr^ 
* dt' 



exp 



0' 



i^Vt - 1' 



r2 



exp 



A{t-t') 

{n-i? 
4{t-t') 



d^ dr 

di + 0{£^') 



In the left-hand side, we change the variables (r — = z and obtain 



C(ri + £z)uj{z) 



dt' 



- 1' 



r2 



exp 



ri 



4{t-t') 



dzd^. 



It is clear that to prove relation (67), it suffices to prove that, with accuracy 
up to small values, the term ez can be omitted in the exponent. 
We consider the expression 

"* dt' 



J 



C{ri + £z)uj{z) 
exp 



Ht - 1' 

{ri-i + ezf 



exp 



dzd^. 



4{t -t') J ^ V 4(t - 1') 

Since the function uj decreases fast, we can assume that \z\ < 6 G (0, 1). 
Next it suffices to consider the integral over t' from to t — e^"*^, because if 
t' € [t — e^~^,t], then the obtained integral can be estimated as 0{e^^~^^^'^). 

In the remaining integrals, we can use the same method as in the proof of 
Lemma m Namely, we transform the difference of the exponential functions: 

{ri-C + £z)^\ [ {ri-if^ 



exp 



4(t-t') 



exp 



4(t-t') 



exp 
7 G (0,1/2), 



(2(r, - - ez)£z 
A{t-t') 



1 



ez 
t-t' 



ez 
t-t' 
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and apply the estimate 

2(r. - F\ - fz\f7.\ \ //t-t'\''' 

< 1, 



V V 4(f-f) ) ) I \ ez 

which, obviously, is valid in the required range of variables. Hence, as in 
Lemma 4, we obtain the estimate 

J = 0(e^), /iG (0,1/2), 

which proves relation ([7T]) . □ 

Similarly to Lemmas H] and El we prove the following statement. 

Lemma 6. The following relation holds for any of the functions qi defined by 
relations 13 1\}: 



Un,t) - qi{m,t) = 0{\i>oH^), /iS (0,1/2). 

Lemma 7. Relations 177] ), [J^ ) hold. 

Proof. In Section 3, we introduced the function 0,{z,rj). Of course, we can 
change the signs in the arguments using the fact that the function woiz) is 
odd, but the form presented above is convenient for calculations. 
Calculating the derivatives ut and Ux, we obtain 

- [/3(ri - r)\ coo (/?^) ^0 (72) 



...|{_.„(,ri^),.„(,^) p3, 



Now we group the terms in the subintegral expression so that one group 
contain terms with the factor uJo(f3{ri — r)/e), and the other group contain 
terms with the factor uJo{l3{r — r2)/e). Thus, we write the integral as the sum 
of two integrals. Then we perform the change of variables 

r^n-j (74) 
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in the first integral and the change of variables 



ez 

To H 

/3 



(75) 



in the second integral. 

From (j72p and (j73p . we obtain 



0{e 



r2t + -^fo 



(76) 



C(ri)[(l + zm -n'^) + {l-z- vKW^ - dz 



C(r2)[(l + z)f^L + r\l - z - - n^] dz 



Here we used the formula for the derivative /?j = Pri^'ot ^^^^ 
that the functions in the integrand in (j76p are invariant under the change 
z — > —z — 7], since each bracket is multiplied by —1 in this change. 

Now we calculate the other expressions contained in Q have 



rX {r){ur) dr 



8e 



r — r2 

e 



n — r 



(77) 



^(rfC'(n)+riC'(r2)) / {d:,{z){l-u^,{-z-ri)f 



- loq{z)loq{-z - r/)(l - ujq{z)){1 - uq{-z - 7]))] dz + 0{e). 

It is easy to see that the integrand in the right-hand side of (j77p can be written 
as 



{{loq{z){1 - loq{-z - 



d;o(z)tj(-z -?/)(!- wo(z))(l - ujq{-z - rf))} dz = ^ j {n'^f dz. 
Thus, we finally obtain 



rWurfdx = ^(rlC'{n)+4C'{r2)) I {Q'^f dz + 0{e). 
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Then the left-hand side of (j77p is positive, and hence the integrand expression 
in the right-hand side is also positive (until the measure of the points at which 
\ur\ ^ const is positive). 
Similarly we obtain: 

i I r^F{u)Cdr = ^{rlC'in) +rlar2)) j F{n)dz + 0{e). 

To prove this relation, we must successively perform the change of variables 
(j75p in the integral in the left-hand side and then consider the half-sum of the 
integrals obtained, which, obviously, is equal to the original integral in the 
left-hand side. Moreover, we must take into account the estimates (j37p . (j38p . 
and 

'l C'(exp(-2z)), z ^ oo, 

uJo{—z — ri) + 0{exp{2z)), z ^ —oo, 

which together with the explicit form of the function F{u) = ^ — ^ + i 
imply that the integrals 



n{z,ri) 



zF{n) dz 

converge. □ 
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